Determination and Measurement of Melanopic Equivalent Daylight (D65) Illuminance (mEDI) in the Context of Smart and Integrative Lighting

In the context of intelligent and integrative lighting, in addition to the need for color quality and brightness, the non-visual effect is essential. This refers to the retinal ganglion cells (ipRGCs) and their function, which were first proposed in 1927. The melanopsin action spectrum has been published in CIE S 026/E: 2018 with the corresponding melanopic equivalent daylight (D65) illuminance (mEDI), melanopic daylight (D65) efficacy ratio (mDER), and four other parameters. Due to the importance of mEDI and mDER, this work synthesizes a simple computational model of mDER as the main research objective, based on a database of 4214 practical spectral power distributions (SPDs) of daylight, conventional, LED, and mixed light sources. In addition to the high correlation coefficient R2 of 0.96795 and the 97% confidence offset of 0.0067802, the feasibility of the mDER model in intelligent and integrated lighting applications has been extensively tested and validated. The uncertainty between the mEDI calculated directly from the spectra and that obtained by processing the RGB sensor and applying the mDER model reached ± 3.3% after matrix transformation and illuminance processing combined with the successful mDER calculation model. This result opens the potential for low-cost RGB sensors for applications in intelligent and integrative lighting systems to optimize and compensate for the non-visual effective parameter mEDI using daylight and artificial light in indoor spaces. The goal of the research on RGB sensors and the corresponding processing method are also presented and their feasibility is methodically demonstrated. A comprehensive investigation with a huge amount of color sensor sensitivities is necessary in a future work of other research.


Introduction
The development of lighting technology from the beginning of the 20th century to the present can be roughly divided into three stages, some of which overlap. Until the end of the 20th century, with thermal radiators and discharge lamps as light sources, the first component of integrative HCL ("Human Centric Lighting") lighting technology was about "visual performance" (e.g., reaction time, contrast perception, reading speed, visual acuity) to enable visibility, improve work performance, and minimize error rates in industrial, educational, and healthcare facilities and offices, among others [1][2][3][4][5][6][7][8][9]. Since the 1970s and even more so since the mid-1990s, with the growing importance of education and the information society, the psychological effects of light (e.g., scene preference, room preference, user satisfaction, spatial perception, attractiveness, and color preference) have become more and more important. As a result of this research process, international and national standards (e.g., EN 12464 standard [9]) have been formulated for lighting design In recent years, there has been a large amount of research and publications on nonvisual lighting effects that have been analyzed and reviewed in books and literature reviews [20][21][22][23]. In the majority of scientific publications cited in these reviews, the input parameters describing lighting conditions were photometric metrics (illuminance, luminance) and/or color temperature. In addition to the optimal lighting design of buildings and the development of modern luminaires according to the criteria of visual performance and psychological lighting effects, it is necessary to define suitable input parameters for the description of non-visual lighting effects to determine their optimal values and to measure these parameters in a wide variety of applications in the laboratory and in the field. To achieve the above goal, the following three research questions for lighting research must be formulated in general terms:

1.
Which input parameters can be used to describe the non-visual effects of light in their variety of manifestations (alertness, sleep quality, hormone production and suppression, and phase shift)? 2.
What values of these input variables are currently considered in the literature to be minimal, maximal, or optimal? 3.
Which measuring devices, sensor systems, and measuring methods can be used to measure the input quantities for non-visual lighting effects and process them in the context of smart lighting in the course of the control and regulation of LED or OLED luminaires on the basis of the definition of personal and room-specific applications?
To define the parameters of non-visual lighting effects (Research Question 1), we must consider that, in the last 20 years, with more activities in the time range between 2016 and 2022, numerous international scientific discussions and analysis tasks have been carried out to specify the right metric for non-visual effects using datasets from experiments of different research groups on nocturnal melatonin suppression as a validation basis. Rea et al. defined and improved the CS model (Version 2018) in the two years 2020 and 2021, wherein they took into account the exposure time t (in hours) and the visual field, as well as modeled the contributions of the ipRGC channel, the S cones, the rods, and the (L + M) channel [24][25][26]. This improved formulation was validated using melatonin suppression datasets [26]. The measurement method for the CS parameter in the 2018 version was published by Truong et al. [27]. As well, the International Commission on Illumination (CIE) [28] and scientists in the fields of neurophysiology, sleep research, and lighting technology have recently made efforts to find effective parameters based on the evaluation of the five photoreceptor signals (LMS cones, rods, and ipRGC cells) and the corresponding calculation tool [29]. This scientific process resulted in the definitions of "Melanopic Equivalent Daylight (D65) Illuminance (mEDI)" and the "Melanopic Daylight (D65) Efficacy Ratio (mDER)". These definitions have been recognized and used by international experts for several years and form the basis of Section 2 of this paper. Additionally, similar issues have been addressed in the WELL Building Standard v2, Q3 2020 version [30,31], which includes recommendations for melanopic equivalent daylight illuminance (mEDI ≥ 136, 109, 218, or 180 depending on the type of light sources). To answer Research Question 2, we recognize that international sleep researchers and neuroscientists made some recommendations in 2022 regarding non-visual light effects, physiological aspects, sleep, and wakefulness based on their literature reviews [32]. The physiological aspects included hormonal regulation cycles, heart rate, core body temperature, and certain brain activities. The main findings were as follows: (a) "Daytime light recommendations for indoor environments: Throughout the daytime, the recommended minimum melanopic EDI is 250 lux at the eye measured in the vertical plane at approximately 1.2 m height (i.e., the vertical illuminance at eye level when seated)". (b) "Evening light recommendations for residential and other indoor environments: During the evening, starting at least 3 h before bedtime, the recommended maximum melanopic EDI is 10 lux measured at the eye in the vertical plane approximately 1.2 m height. To help achieve this, where possible, white light should have a spectrum depleted in short wavelengths close to the peak of the melanopic action spectrum". (c) "Nighttime light recommendations for the sleep environment: The sleep environment should be as dark as possible. The recommended maximum ambient melanopic EDI is 1 lux measured at the eye".
The first recommendation, with a minimum value of "melanopic EDI = Melanopic Equivalent Daylight (D65) Illuminance, mEDI" of 250 lx measured vertically at the observer's eye, is relevant and of great interest for the professional sector during the day (offices, industrial halls, educational and health facilities). The other two recommendations are for the dark hours in the home. For the solution to Research Question 3 "Determination and measurement of input parameters for non-visual lighting effects", which builds the focus of this present paper, two application areas can be targeted: 1.
The determination and measurement of non-visual parameters after the completion of new lighting installations and comparison with the specifications of the previous lighting design; or verification of the results of the development of new luminaires for HCL lighting in the lighting laboratory or in the field. For this purpose, absolute spectroradiometers for spectral radiance or spectral irradiance are used to calculate the parameters mEDI and mDER at different locations in the building. These two non-visual parameters cannot be measured directly with integral colorimeters and illuminance-luminance meters.

2.
The control or regulation of modern semiconductor-based lights (LED-OLED) and window systems (daylight systems) with the help of sensors in the room (e.g., presence sensors, position sensors, light and color sensors) [33]. In order to achieve a predefined value of non-visual parameters such as mEDI at a specific location in the room (e.g., in the center of the room or at locations further away from the windows), while taking into account dynamically changing weather conditions and the whereabouts of the room users, in practice, relatively inexpensive but sufficiently accurate optical sensors (RGB sensors) are required. The goal is to obtain not only the target photometric and colorimetric parameters such as illuminance E v (in lx), chromaticity coordinates (x, y, z), or correlated color temperature (CCT in K), but also the non-visual parameters mEDI and mDER. The principle of this Smart Lighting concept using RGB color sensors is illustrated in Figure 2. The methods for measuring non-visual parameters with low-cost but well-qualified RGB color sensors are the content of the Sections 3 and 4 and the focus of this paper. The notations "integrative lighting" and "human-centric-lighting" used in this publication have been defined by the ISO/CIE publication [34] and mean a lighting concept and practice that integrates both visual and non-visual effects and produces physiological as well as psychological benefits for light users. The terms "smart lighting" or "intelligent lighting" are used in the same way and mean a framework that combines the above-mentioned integrative lighting as a goal and content of lighting practice with the technological aspects of lighting (signal communication, Internet of Things, sensor systems, control electronics, software that includes the methods of artificial intelligence, LED luminaires) while taking into account the individual needs of light users. In order to achieve the above objectives, the content of this paper is structured as follows: • Section 2 defines the non-visual parameters according to the CIE publication [28].
• Section 3 describes the characterization and signal transformation mathematics for converting the RGB signals of the RGB sensors to the chromaticities x, y, and z according to the CIE publication [35] for a viewing angle of 2°. • In Section 4, a computational model is established and analysed that converts the chromaticity z into the non-visual parameter mDER and the combination of z and illuminance E v into mEDI. Section 5 summarizes the results with some outlook.

Definition of Non-Visual Input Parameters [28]
The extensive mathematical treatment and definition of the non-visual input parameters are described in detail in the CIE publication S 026/E:2018 "CIE System for Metrology of Optical Radiation for ipRGC-Influenced Responses to Light" [28], and those contents are summarized in this paper. The purpose of the CIE publication is to define spectral response curves, quantities, and metrics to describe optical radiation for each of the five photoreceptors in the human eye that may contribute to non-visual processes. The spectral sensitivity functions of the five receptor types are specified as follows (see Figure 3): 1.
Equation (5) is similar to Equation (3.7) on Page 5 of [28]. For melanopic efficacy (α = mel), the parameter value K D65 mel,V equals 1.3262 mW/lm. Consequently, the parameters mEDI and mDER can be set up as shown in Table 3. If one wants to interpret the meaning of mEDI or mDER for lighting engineering, there are two main aspects: mEDI is the illuminance of the standard daylight illuminant D65 that has as much melanopic efficacy as the test light source for a given illuminance E v (lx) caused by the test light source. See Equation (6) in Table 3.

2.
mDER is the ratio of the illuminance of the standard illuminant D65 (mEDI) to the illuminance of the test illuminant E v (in lx) when the absolute melanopic efficacy of both illuminants is set equal. See Equation (9) in Table 3.
The melanopically effective parameters mDER and mEDI can be calculated from the Equations (3), (5), (6), and (9) if the spectral irradiance E e,λ (λ) is known from a spectroradiometric measurement. The key question in this publication is whether these two parameters can also be determined with sufficient accuracy in practical lighting technology using a well-characterized and inexpensive RGB sensor in the sense of intelligent lighting technology (smart lighting). Section 3 deals with the description of the RGB sensor and how the RGB sensor signals can be transformed into tristimulus values (RGB) and into chromaticity coordinates (x, y, z) by means of a comprehensive spectral analysis.
(*) Note: For the parameter K y x (indices x, y according to the corresponding definitions), see Equations (4) and (5) when α = mel, and apply Equation (3) when the calculated parameter is the mel-opic irradiance.

RGB Color Sensors: Characterization and Signal Transformation
After the definition of the non-visual input parameters [28] in Section 2, we can understand the main concepts and mathematical forms of mEDI and mEDR. In Section 4, the prediction of a simple computational model for the non-visual quantities mEDI and mDER was carried out. As well, the verification of the feasibility of a RGB sensor using this model was implemented to check the synthesized prediction model for lighting applications. This was also performed to confirm that the model not only conformed to the mathematical calculations, but also to determine whether it can be applied in the lighting systems with low-cost RGB color sensors. Therefore, Section 3 must attempt to describe RGB color sensors, their characterization, and appropriate signal transformation techniques. This section serves as a link between Sections 2 and 4, as well as to balance the essential material for the verification of the prediction model using a color sensor as an example in Section 4 later. As a demonstration of the methodology for processing RGB sensors, it is not necessary to collect all color sensor sensitivities, synthesize, and compare them, but the most important thing here is to prove that the methodology can work well in this approach. Consequently, the future work can be implemented more comprehensively for different color sensor families through other research.

Characterisation of RGB Color Sensors
The color sensors (see Figure 4) include arrays of individual sensors or a group of sensors based on semiconductors with silicon as the semiconductor material, which are most commonly used for the visible spectral range between 380 nm and 780 nm. In order to generate the corresponding spectral sensitivity curves of the individual color channels in the red, green, and blue range (so-called RGB channels), optical color filters based on thin-film technology (interference filters), absorption glass filters or color varnishes were microstructured and applied to the respective silicon sensor. For the correct spatial evaluation of the optical radiation according to the cosine law, a so-called cos prefix (usually a small diffuser plate made of optical scattering materials) was applied to the RGB sensor. The photons of each wavelength are absorbed by the RGB sensor, and photon currents are generated, which are then converted into voltages in an amplifier circuit. These voltages are then digitized in different bit depths (8 bits, 12 bits, 16 bits) by an A/D converter (analog-digital). The spectral sensitivity of each color channel R-G-B is, therefore, made up of the components shown in the Equations (10)- (12). where: • S(silicon, λ) is the spectral sensitivity of the silicon sensor; • τ(di f f usor, λ) is the spectral transmittance of the cos diffusor; • τ(color lacquer i , λ) is the spectral transmittance of the color filter layer for the respective color channel i = R, G, B; • K R , K G , and K B are absolute factors to account for current-to-voltage conversion and voltage digitization; • CV i (λ) values are output signals (analog or digital) of the respective R-G-B color sensor at wavelength λ.
If a certain spectral irradiance E(λ) is present on the RGB sensor, three corresponding output signals R, G, and B are generated in the three color channels R, G, and B, respectively (see the Equations (13)-(15)).
For the calculation of the RGB channel signals, it is, therefore, necessary to know or to determine the spectral sensitivity function of each individual color channel R, G, B by laboratory measurements. Consequently, the spectral apparatus for determining the spectral sensitivity of the semiconductor sensors in the authors' light laboratory (see Figure 5) consists of a high-intensity xenon ultrahigh-pressure lamp, a grating monochromator (spectral half-width ∆λ = 2 nm, spectral measuring steps for the entire spectrum between 380 nm and 780 nm, ∆λ = 2 nm), and an integrating sphere for homogenizing the quasi-monochromatic radiation coming out of the monochromator. The RGB color sensor and a known calibrated reference sensor are located at two different locations on the inner wall of the sphere behind a shutter. During spectral measurement of the RGB sensor, one can determine the spectral sensitivity for the red color channel, for example, according to Equation (16).
with CV R (λ) as the output signal of the red color channel and E(sphere, λ) as the spectral irradiance at the inner wall of the sphere when the monochromator is set to the wavelength λ, which in turn can be determined from the measured output photocurrent i(reference, λ) and the known spectral sensitivity of the reference sensor S(reference, λ). As a result, we can write Equation (17).
The absolute spectral sensitivity of each color channel R, G, B can be determined from the Equations (16) and (17). Figures 6 and 7 show examples of the spectral response curves of some RGB color sensors measured in the authors' light laboratory. In Figure 6, the two datasets of the same type, SeS1 and SeS2, differ by the different peak heights of the R channel and by the different slopes of the G and B curves, because the two sensors SeS1 and SeS2 came from different production sets. The SeS1 and SeS2 spectral sensitivity curves in Figure 6 are fundamentally different from the RGB spectral sensitivity curves of the color sensor type in Figure 7.  The most accurate characterization possible for RGB sensors and the best possible similarity between CIE color matching functions CIE xyz(λ) for a field of view of 2°and RGB(λ) sensitivities are very important and have a direct influence on the later determination of non-visual parameters such as the results in Table 13. The smaller the difference between RGB(λ) sensitivities and CIE color matching functions CIE xyz(λ) for a field of view of 2°, the closer the distance between the values RGB SeS. , the chromaticity coordinates xyz CIE , and later the values xyz SeS. after matrix transformation of the chromaticity coordinates xyz CIE . Finally, the difference between the determined z SeS. = (1 − x SeS. − y SeS. ) and the CIE chromaticity coordinates z CIE. = (1 − x CIE. − y CIE. ) must be as small as possible, and this is understood as the primary uncertainty of the non-visual parameter calculation. Then the non-visual parameter model mEDR = f (z SeS. ) in Section 4 and the parameter mEDI = [mEDR(z SeS. ) · E v ] determine the secondary uncertainty of the non-visual parameter calculation. The abbreviation "SeS." indicates that it is the value of the sensors. Figures 6 and 7 show that the RGB curves of the real color sensors differ more or less strongly from the xyz curves of the CIE color matching function for a field of view of 2° [35]. In order to obtain the XYZ tristimulus values, the generated RGB signals must be transformed into digital form by matrixing (Equation (18)).

Method of Signal Transformation from RGB to XYZ
the sensor technology, the amplification electronics, and the A/D conversion often have a more or less pronounced non-linear behavior, it may be necessary to use different matrixing types from 3 × 3 to 3 × 22. See Table 4. Table 4. Different matrix types for color sensing according to [36].

Nr.
Size Content If the sensor electronics are linear and the spectral response curves have a similar relative shape to the xyz color matching functions or LMS sensitivity spectra of the retinal photoreceptors, the 3 × 3 to 3 × 8 matrices will yield contributions from RGB signals as a function of the first term. The more the sensor electronics deviate from linear behavior, the more the shape of the spectral curves of the real sensors deviate from the xyz color matching functions. Therefore, matrices with RGB contributions in quadratic or cubic functions (see Table 4) should be set up. The procedure for finding the optimal matrix based on the chromaticity difference ∆u v is shown in Figure 8.
The goal of the optimization is to get the chromaticity (x, y, z) of the RGB sensor as close as possible to that derived from the CIE calculation. Instead of working with the non-uniform xy color diagram, it is better to work with the more uniform u v diagram (CIE, 1964). Bieske [37] showed in her dissertation that, when ∆u v is less than 10 −3 , the color difference is not perceived by the subjects. From 10 −3 to 3 × 10 −3 , the color difference can be perceived but is still acceptable. If this value is higher than 5 × 10 −3 , the color difference is unacceptable to the subjects. Based on this scientific contribution, the optimization with u v not only gives the best z-value, but can also be directly checked and compared with the obtained perception thresholds.

Matrix Transformation in Practice and Verification with a Real RGB Color Sensor
For the color sensor type whose spectral sensitivity curves are shown in Figure 7, the optimization processes are carried out according to the scheme in Figure 8 with a series of matrix types from 3 × 3 to 3 × 22 with nine different lamp spectra (as a training set). These nine lamp spectra are tabulated in Table 5 and shown in Figure 9. The illuminance for the optimization of the matrix was chosen to be E v = 750 lx. In this training set, the following lamp types were selected to reflect the variety of types used in practice today: fluorescent lamps between 2640 K and 4423 K with many spectral lines (see Figure 9), a halogen incandescent lamp, two phosphor-converted white LEDs, and a combination of RGB LEDs and warm white LEDs (RGBWW4500K). The calculation with different matrix types and with individual lamp spectra yields color differences ∆u v of varying magnitudes, with the 3 × 3 matrix form yielding the smallest color difference (see Table 6). Optimization with more lamp types did not bring any significant improvement in this case. The final 3 × 3 optimized matrix for the transformation from RGB to XYZ, as well as the formula for the calculation of the illuminance from the RGB signals, are shown in Table 7.    To verify the prediction quality of the formula in Table 7, eight spectra were selected in the next step. These were white phosphor converted LEDs (pc-LEDs) and a mixture of daylight spectra with white phosphor converted LEDs (TL-pc-LEDs), as is often the case in practical indoor lighting (see Figure 10). The illuminance was again set to 750 lx. The results of the check are shown in Table 8. From Table 8 it can be seen that The deviation of the illuminances, calculated directly from the lamp spectra via the RGB color sensor signals and via the formula in Table 7, was below 0.65% in percentage terms; (b) For the majority of pc-LEDs and for the combinations daylight-white LEDs, the chromaticity difference ∆u v was in a small or moderately small range from the point of view of practical lighting technology, although the matrix transformation was synthesized only for the general case of many lamp types. No specific procedure was optimized for the specific case of only mixed light between daylight and conventional light, but the color and illuminance difference was still very small when using the achieved matrix transformation. Therefore, the authors did not build a separate case of mixed light with specific matrix transformation. Also, pc-LEDs and their mixing with daylight are dominant nowadays, so this verification in this way makes sense for applications. Therefore, the extreme cases of very complicated spectral shapes are not necessary to investigate their feasibility, because they are very rarely applied in today's life. And an exception is the spectrum R12GB12-5000 K (green curve in Figure 10) with three distinct peaks of the three RGB LEDs (B around 455 nm, G around 525 nm, R around 660 nm), which is a very good demonstration example of a complicated spectral form. In this case, the color difference was ∆u v = 10.1 × 10 −3 . However , for the case of pure daylight spectra, a special matrix was found and described in Table 9 as a specific case. For the investigated RGB color sensor type, 185 daylight spectra were measured absolutely on a sunny summer day (19 August 2020 in Darmstadt, Central Europe) from 6:32:00 in the early morning to 20:35:32 in the late evening (see Figure 11). The chromaticity and illuminance values were calculated and used as a training set for the matrix. The largest color difference with a 3 × 3 matrix was found to be ∆u v = 1.2 × 10 −3 , which was recorded in the last minutes of the evening before sunset when the correlated color temperature was very high, measuring approximately in the range of 17,000 K on the day of measurement. Most of the color differences were well below this value.  Table 7. Table 8. Verification of optimal matrix (3 × 3) in Table 7 with 8 LEDs plus mixed daylight spectra.
To verify the 3 × 3 matrix for daylight, 24 daylight spectra measured on a overcast day (23 September 2020 in Darmstadt, Germany) were considered. The verification results are shown in Table 10 with 9 out of 24 spectra as examples. There, the chromaticities and illuminances determined from the measured spectra are compared with the data obtained by matrixing the RGB color sensor (processed data). The calculated color difference ∆u v was very small and lay in the range of 10 −4 . The determination of the optimal matrix for the transformation of RGB sensor signals into XYZ values and the illuminance E v described so far, as well as the verification results with an actual RGB color sensor type, prove that it is possible to obtain the colorimetric and photometric values XYZ, E v , and CCT with relatively good results in the sense of a reliable, relatively accurate and adaptive lighting technology using inexpensive and commercially available color sensors. The next section deals with the accurate determination of the nonvisual input quantities mEDI and mDER indoors (daylight and artificial light combined or separated) and outdoors during the day (with daylight only) from the tristimulus values XYZ and the illuminance E v of a qualified RGB color sensor, because, in practical lighting technology, a spectroradiometer is often too expensive and too cumbersome to handle.

Prediction of a Simple Computational Model for the Non-Visual Quantity mDER and Verification of the Feasibility of a RGB Sensor Using this Model
Once the matrices for electric light sources and for daylight spectra have been found, the tristimulus values XYZ, the chromaticity coordinates xyz, the color temperature CCT, and the illuminance E v can be obtained from the RGB color sensors used (see the  Tables 7-10). To obtain the melanopic equivalent daylight illuminance mEDI, one only needs to know the value of mDER and the illuminance E v (in lx) according to Equation (7), where the illuminance can be formed from the RGB signals (see the formula in Tables 7 and 9) for determination. Consequently, the most important task is to calculate the dimensionless input quantity mDER (see Equation (9)) from the chromaticity coordinates x, y, z (z turned out to be the most predictive, so it was the coordinate used).
For this transformation from z to mDER, a large number of measured and simulated light source spectra were analyzed (see Table 11). The first four light source groups were real measured light sources from thermal radiators (28 light sources), compact and linear fluorescent lamps (252 light sources), different LED configurations (419 light sources), and 185 measured daylight spectra. These 884 light source spectra are shown in Figure 12. In order to simulate the combination of daylight with white LEDs or daylight with white fluorescent lamps, which can correspond to the lighting conditions of interiors with white LEDs or white fluorescent lamps with windows over a longer period of use, the authors of this paper took 185 measured daylight spectra and mixed them with a market-typical white LED spectrum or with a typical fluorescent lamp spectrum in nine different mixing ratios from 10% to 90%. This can be seen in Figure 13 for the case of combining daylight spectra with white LEDs. For the fluorescent lamps, it was similar, so a graphical representation of the mixing spectra is omitted here. The 185 daylight spectra with nine mixing ratios each resulted in 1665 simulated spectra each for the LED and fluorescent lamps (185 × 9 = 1665).   Using the 4214 light source spectra summarized in Table 11, 4214 mDER values and 4214 z values were calculated, from which their relation was determined. The formula for mDER is given in Equation (25).
The parameters a, b, c, and d for the fitting function and the correlation coefficients R 2 are given in Table 12. The course of the correlation between mDER (ordinate) and the chromaticity coordinate z (abscissa) is shown in Figure 14, where the fitting function follows the green course of the curve.

Figure 14.
Correlation between mDER (ordinate) and the chromaticity coordinate z (abscissa). To check the predictive quality of Equation (25) and the transformation from z to mDER, the 4214 spectra in Table 11 were first used to calculate the 4214 mDER values from the spectra themselves (denoted as mEDR [CIES026/E:2018] ) and then another 4214 mDER values indirectly via Equation (25), with z as the independent variable (denoted as mEDR = f (z)) (see Figure 15). The goodness of fit (R 2 ) was 0.97844, and the linear constants were 1.0214 and −0.05712. Compared to the ideal constants of 1 and 0, the quality of the fit is very good. To verify the quality of the characterization of the RGB color sensor under test (see Figure 7), the quality of the "matrixing" in the transformation from RGB to XYZ and E v (see Table 7) and the quality of the transformation from z to mDER (see Equation (25)) and finally the calculation of mEDI (see the Equations (6) and (9)), we set each of the eight light source spectra in Figure 10 to 750 lx. From these 8 spectra, the chromaticity coordinates x and y and the quantities CCT, CRI R a , E v , mDER original , and mEDI original were calculated. They are listed in Table 13. For these eight spectra, the chromaticity coordinates x new and y new , the value of E v,new , and the values mDER new and mEDI new were also calculated from the RGB color sensor values by matrixing and transforming z to mDER and mEDI, respectively. The maximum relative deviations ∆mDER (in %) and ∆mEDI (in %) were found to be in the range of ± 3.3%.

Conclusions and Discussion
In indoor lighting technology, including daylight components during the day, the daily task is to evaluate the lighting systems after the completion of the building or after reconstruction to assess the photometric and colorimetric quality of the newly developed interior luminaires and, in the context of intelligent lighting (smart lighting, HCL lighting, integrative lighting), to adaptively control the lighting systems in order to provide the room users with the best visual and non-visual conditions, as well as room atmosphere, at all times. In addition to the criteria of visual performance according to current national and international standards ( e.g., EN DIN 12464 [9] or WELL [31]), aspects of psychological and emotional lighting effects, as well as the evaluation and adaptive control of lighting systems according to non-visual quality characteristics, should be considered. For this purpose, non-visual input parameters such as mEDI and mDER need to be measured and processed with sufficient accuracy, reliability, and reasonable effort. Existing illuminance meters, luminance meters, color and luminance cameras, and small portable colorimeters are currently only capable of measuring photometric quantities such as illuminance E v , luminance L v , colorimetric parameters such as chromaticity coordinates xyz, and correlated color temperature (CCT). To measure the non-visual parameters, a transformation is needed to convert the tristimulus values to the non-visual parameters, and a sensor platform is needed to convert the RGB sensor signals to the tristimulus values XYZ and E v (in lx). The authors of this paper have characterized some exemplary RGB color sensors on the basis of laboratory measurements, calculations, and optimizations by creating one or more matrices for the optimal transformation of RGB sensor signals to the tristimulus values. For different spectral sensitivities of the RGB color sensors, different matrices with linear or higher order RGB-signals can be optimized. As an example, for a certain RGB sensor, a 3 × 3-matrix was found by means of optimization with nine different spectra to be the best transformation that delivered a maximal chromaticity difference of ∆u v = 0.0083 (see Table 6). This matrix was verified with eight white phosphor-converted LEDs (pc-LEDs) and mixture of daylight spectra with white phosphor-converted LEDs (TL-pcLEDs) (see Table 8). For the spectra of daylight only, a special matrix was found that allowed for a much lower chromaticity difference (maximal ∆u v = 0.000857, see Table 10). The reconstruction of the illuminance (in lux) from the RGB signals delivered a difference in the order of 0.64% (see Table 8). As next steps, on the basis of extensive calculations, a transformation from the chromaticity coordinate z to mDER and, via E v , also to mEDI, was achieved with good accuracy. Table 13 has shown the verification results with eight light source spectra, which were in the range of ∆mEDI = ±3.3%, which can be accepted from a practical point of view of lighting engineers, lighting designers, and luminaire system developers. This work and the methodology described therein enable an accurate and financially justifiable assessment of lighting installations according to non-visual criteria, which will become increasingly important in the coming decades. This paper followed the aim to deliver, firstly, a framework and a general method on how to use commercially available RGB sensors to qualify them for measuring the metrics for non-visual effects. The results had some limitations, because only one RGB sensor type has been used in the methodical frame of this paper. Other RGB sensor types may have other spectral, optical, and electronic characteristics with different dark currents, signal-to-noise ratios, and cross-talk effects, which should require other forms of matrices (e.g., matrices with higher orders of RGB signals). Further studies will be planned by the authors of this present paper, in which the accuracies of several commercial RGB sensor types will be compared to the measurement results of an absolute measuring spectroradiometer for different test conditions (e.g., outdoor and indoor lighting conditions with different mixture ratios of daylight and artificial light sources). Table 13. Verification of the mDER and mEDI values with the color sensor under test on the basis of the 8 spectra in Figure 10 set to 750 lx Name pcLED-WWG3 pcLED-WWVG1 pcLED-WWEX4 pcLED-NWVG1 pcLED-NWEXCV2 pcLED-CWEX2 R12GB12-5000 K TL-pcLED-30%  Funding: This work received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors. The publication of the manuscript is supported by the Open Access Publishing Fund of the Technical University of Darmstadt.